Effect of accretion on primordial black holes in Brans-Dicke theory 
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We consider the effect of accretion of radiation in the early universe on primordial black holes in 
Brans-Dicke theory. The rate of growth of a primordial black hole due to accretion of radiation in 
Brans-Dicke theory is considerably smaller than the rate of growth of the cosmological horizon, thus 
making available sufficient radiation density for the black hole to accrete causally. We show that 
accretion of radiation by Brans-Dicke black holes overrides the effect of Hawking evaporation during 
the radiation dominated era. The subsequent evaporation of the black holes in later eras is further 
modified due to the variable gravitational "constant", and they could survive up to longer times 
compared to the case of standard cosmology. We estimate the impact of accretion on modification of 
the constraint on their initial mass fraction obtained from the 7-ray background limit from presently 
evaporating primordial black holes. 

PACS numbers: 98.80.Cq, 97.60.Lf, 04.70.Dy 
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I. INTRODUCTION 

The Brans-Dicke (BD) theory [jj which was proposed 
in 1961 is regarded as a viable alternative of Einstein's 
general theory of relativity (GTR). In the BD theory, 
the value of gravitational constant is set by the inverse 
of a time-dependent scalar field which couples to grav- 
ity with a coupling parameter lu. GTR can be recovered 
from this BD theory in the limit of a; ^ 00. The BD 
theory has been used in attempts to understand many 
cosmological phenomena such as inflation early and 
late time behaviour of the universe [1], cosmic acceler- 
ation and structure formation and the coincidence 
problem^. It is also well known that BD-type models 
arise as low energy effective actions of several higher di- 
mensional Kaluza-Klein and string theories Q. 

Black holes which could be formed in the early Uni- 
verse through a variety of mechanisms are known as pri- 
mordial black holes [PBHs]. Some of the well-studied 
mechanisms of PBH formation include those due to infla- 
tion initial inhomogeneities [l,[io|, phase transition 
and critical phenomena in gravitational collapse [ill [T3l . 
bubble collision [l^ or the decay of cosmic loops |14| | 
The formation masses of PBHs could be small enough 
for them to have evaporated completely by the present 
epoch due to Hawking evaporation p^ . Early evaporat- 
ing PBHs could account for baryogenesis 
universe. On the other hand, longer lived PBHs could 
act as seeds for structure formation or even as precursors 
to supermassive black holes observed presentlyfl8|. Fur- 
thermore, in certain scenarios it is possible for the PBHs 
to survive till date and form a significant component of 
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dark matter [13|. It has been shown recently that PBHs 
in the braneworld scenario can efficiently accrete radia- 
tion [13, [m H^l making them considerably long-lived. 

The possibility of black hole solutions in BD theory 
was first proposed by Hawking [23J. Using scalar-tensor 
gravity theories Barrow and Carr [2J] have studied PBH 
evaporation during various eras. It has been recently 
observed that in the context of generalised Brans-Dicke 
theory, inclusion of the effect of accretion leads to the 
prolongation of PBH lifetimes [1^ . The coexistence of 
black holes with a long range scalar field in cosmology 
could have interesting consequences [2^. The possibility 
of variation of fundamental constants of nature over cos- 
mological scales is a fascinating albeit contentitious issue 
, and the black holes themselves could be used to con- 
strain the variation of fundamental constants [2^. An- 
other interesting issue of gravitational memory of black 
holes in BD theory has also been studied [29^, . 

Accretion of radiation by PBHs in the radiation domi- 
nated era of the early universe has been a much debated 
issue in standard cosmology. It is widely held that accre- 
tion is ineffective for sufficiently increasing the mass of a 
PBH jsij , though later works have poi nted out contrary 
possibilities in certain cases [13, l32l. |33| . Recently, it has 
been realized in the context of the braneworld scenario 
that the possibility of enhanced accretion is quite favored 
due to the modified PBH geometry, as well as the mod- 
ified early high energy era of the universe [13, [llj . Such 
a feature of effective early accretion prolonging the PBH 
lifetime by significant orders could also be valid for other 
modified gravity scenarios, as has already been shown in 
the context of a generalized scalar-tensor model [25*1 . 

The motivation for the present work is to study how ac- 
cretion of radiation during the radiation dominated era 
impacts the evolution of primordial black holes in BD 
theory. Cosmological solutions during different eras in 
the BD theory were obtained by Barrow 0|. Here we 
consider together the processes of accretion of radiation 
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and Hawking evaporation for PBHs in BD theory using 
the solutions for the scale factor and the gravitational 
"constant" as used by Barrow and Carr In the 

present work we do not consider the effect of gravitational 
memory [2^ Is^] on PBH evolution, but assume that the 
evolution for the BD field is similar for the PBHs as it is 
for the whole universe. Cosmological observations could 
be used to impose constraints on the density of primordial 
black holes at various eras [1^. Within this formalism 
we also estimate the impact of accretion in modifying the 
constraint on their initial mass fraction in BD theory ob- 
tained from the 7-ray background limit from presently 
evaporating primordial black holes. 



II. PBHS IN BRANS-DICKE THEORY 

For a spatially flat(fc — 0) FRW universe with scale 
factor a , the Einstein equations and the equation of mo- 
tion for the JBD field $ take the form 
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we consider accretion which is effective in radiation domi- 
nated era, then the primordial black holes take more time 
to evaporate. Let us study how accretion changes the life 
time of the primodial black holes. For a primordial black 
hole immersed in the radiation field, the accretion of ra- 
diation leads to the increase of its mass with the rate 
given by 



Mace = AnfR^HPR 



(7) 



where Rbh ~ 2MG is the black hole radius, — 
8§G (a) radiation energy density surrounding the 

black hole, and / '--^ is the accretion efficiency. The 
value of the efhciency of accretion / depends upon com- 
plex physical processes such as the mean free paths of 
the particles comprising the radiation surrounding the 
PBHs. Any peculiar velocity of the PBH with respect to 
the cosmic frame could increase the value of / [13] ■ Since 
the precise value of / is unknown, it is customary [2l| to 
take the accretion rate to be proportional to the product 
of the surface area of the PBH and the energy density of 
radiation with / ~ 0(1). After substituting the above 
expressions for Rbh and pu equation ([7]) becomes 



(8) 
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The energy conservation equation is given by 

p + 3{j+l)Hp = (4) 

assuming a perfect fluid equation of state p — jp . 

Barrow and Carr [g^l have obtained the following so- 
lutions for a and G for different eras, as 



t(l-^)/3 {t < ti) 

ait) (X { {ti <t < te) 



(5) 



and 
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where ti is the time of starting of radiation dominated 
era, te ~ is the era of matter-radiation equality, to ^ is 
the present time, Go ^ is the present value of G ~ 4?^, 

and n is a parameter related to oj, i.e., n = . 
Since solar system observations [s^ require that uj be 
large {uj > 10''), n is very smaU (n < 0.00007) . 

Barrow and Carr have considered only evaporation of 
the primodial black holes due to Hawking radiation. If 



Accretion is effective only in the radiation dominated 
era. So the primordial black holes which exist during the 
radiation dominated era are affected by accretion. Using 
the solutions for the scale factor a{t) ^ and G{t) ^ in 
equation ([8]), we get for the radiation dominated era 
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where Mi is the black hole mass at it's formation time 
ti > ti. Assuming the standard mechanism for PBH 
formation due to gravitational collapse of density per- 
turbations at the cosmological horizon scale 3, we have 
Mi ~ G~^{ti)ti . Since the horizon mass grows as 
Mnit) ^ G^^t, one finds that Mh grows faster than 
the black hole mass Mbh which for large times asymp- 
totes to Mi{l — 3/2/)"-'. Thus enough radiation density 
is available within the cosmological horizon for a PBH 
to accrete causally, making accretion effective in this sce- 
nario. However, the maximum accretion efficiency can- 
not exceed / = 2/3, thereby making it improbable for the 
overall mass of a PBH to increase much due to accretion. 
Nonetheless, the occurrence of accretion prolongs the on- 
set of the evaporating era for the PBH thereby prolonging 
its lifetime considerably which in turn could significantly 
impact the observational constraints on PBHs in different 
eras. 

Due to Hawking evaporation, the rate at which the 
primordial black hole mass decreases is given by 



Mevap — —'^'^Rbh'^hTbh 



(10) 
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where an is the Stefan-Boltzmann constant, and Tbh 
g^Qj^j is the Hawking temperature. Now 
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If we consider both evaporation and accretion simulta- 
neously, then the rate at which the primordial black hole 
mass changes is given by 



Mbh = 6.fG 



1 
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This equation can not be solved exactly. But we can very 
well approximate it during different regimes when either 
accretion or evaporation is the dominant process. Subse- 
quently, we also integrate it by using numerical methods, 
which corroborates our approximation. 



III. PBH DYNAMICS IN DIFFERENT ERAS 

A. For t <ti : 
Hawking evaporation rate for this era is given by 



\t^) \TJ VlkP 



(13) 



where a 
gets 



Integrating the above equation, one 



M3 = Mf + 3«(l + 2 V^) (|) (ti' 

(i2v^+l _ i2v^+l^ (14) 

This regime corresponds to the BD field dominated dy- 
namics, where the radiation density is only subdominant. 
Assuming accretion is not effective in this era, we get 
same result as that of Barrow and Carr [23| for the evap- 
oration time (A/ = 0) . 
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B. For ti<t<U 



This period corresponds to the radiation dominated 
era. Here we consider two possibilities: PBHs created 
before ti and those created after <i. 

CASE-I (t, < ti) : 

The PBH evaporation equation (fTT|) becomes 
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Integrating this equation, one gets 

Vto/ 



U " t 



(17) 



Now considering both evaporation and accretion, we 
obtain 

"-—(i^) % (V) 



during the BD field dominated era (t <tx} and 



M.„^|/0.(i^)"(^)-«(|)"'l^ 
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for t» t\ in the radiation dominated era. Since accre- 
tion is effective only during the radiation dominated era, 
Eq. p^ is first integrated over the time period ti to t-y 
and then Eq. (fTO)) is then integrated over the time period 
t\ to t with the initial condition Mbh^^x) obtained from 
the solution of Eq. (fT5|) . The results of numerical integra- 
tion of the above equation with several values of Mi are 
presented in the figure- 1. We assume that radiation dom- 
ination sets in at the GUT scale, i.e., t\ = lO"^^^. We 
find that the PBHs formed before the onset of radiation 
domination, i.e., for ti < ti, evaporate out completely 
during the radiation dominated era. Figure 1 shows the 
variation of PBH masses for three different values of Mi 
with Mi ~ 990(7 corresponding to the formation time 
t ^ ti. Though here the accretion efficiency / = 1/3 , 
it can be verified that even for larger /, accretion in this 
case is not able to prolong PBH lifetimes beyond the ra- 
diation dominated era, i.e., those PBHs which form at a 
time ti < ti, have evaporations times t^vap < ^e- Note 
that though in our analysis we assume radiation domi- 
nation at the GUT scale, it is possible for the radiation 
domination to set in much later, i.e., ti >> lO^^s^ jn 
certain scenarios, e.g., gravitino production could con- 
strain the reheating temperature to be much lower [STj . 
In such a case of low reheating temperature, accretion of 
radiation will be effective for a shorter window of time 
corresponding to a shortened span of radiation domina- 
tion, and the evolution of PBHs would be closer to the 
case in BD theory without accretion [23 |. 

CASE-II (t, > ti) : 

Taking both accretion and evaporation into account, 
we can write 



ft,\2n 1 



M2' 



(20) 



For PBHs with formation mass Mf > °gg'^^ , the mag- 
nitude of the first term (accretion) exceeds that of the 
second term (evaporation). In the radiation dominated 
era for a PBH whose formation mass satisfies the above 
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Accretion efficiency=l/3 




-44 -42 
log(t/ro) 

FIG. 1: Evolution of PBH masses witli time (to being the 
present age of the universe) for different initial mass Mi — 
0.5(7, lOOg, 990(7, but with same accretion efficiency / = 1/3. 



relation, accretion is dominant upto a value of t, say tc at 
which accretion equals evaporation (the PBH mass rises 
to a maximum value Mmax at this stage) , and after that 
evaporation dominates over accretion. For t = tc, the 
magnitude of the accretion term is equal to the magni- 
tude of evaporation term. So for the radiation dominated 
era, equation (|^ implies, 



which gives 



2 

max 



ytj \Ml. 



where A ^ ^Gn^a(^^^ 



3n 



and Mr, 



3--'0 ^\ta J 

from the PBH accretion equation 
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M,. 
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we have 

-1 



(21) 



(22) 



But 



(23) 



Equating above two expressions for M„ 
7 \ 3 M, 



one gets 



t 



1/2 _ 



A 



and 



M„ 



1 



f 



(24) 



1 



/ 



(25) 



which again stipulates that / < |. 

Considering evaporation from tc onwards, we get 



M = M„ 



tc 



So the evaporation time for these PBHs are given by 

I / ''() \ / ^*^-^m( 



te 



tc 



1 



C. For t>te: 

The PBHs which are formed before radiation domina- 
tion completely evaporate out during the radiation dom- 
inated era. So for t > t,,, only those PBHs which are 
formed after ti exist. The PBH evaporation equation 
(fTTj) can be written as 



M 



dM = 



Mi 



(I) 



2n 



dt 



{toy^'H^'^dt 



(28) 

Taking both accretion and evaporation into account, 
one gets 



M 



BH 



during the radiation dominated era and 



BH 



2ri / +2ri 



(29) 



(30) 



during the matter dominated era. In order to obtain the 
PBH lifetime, one can numerically integrate the above 
equations (|^ and (|5D|) . Eq. (|^ corresponding to ac- 
cretion and evaporation during the radiation dominated 
era is integrated over the period ti to tc and the value of 
MBHite) is obtained and used as the initial condition for 
integrating Eq. ((30|) over the period tc to t. The results of 
numerical integration are displayed for a particular initial 
mass in figure-2. (In the figure we do not show the early 
part of their evolution where their mass increases by a bit 
due to accretion and then stays nearly constant for a long 
period of time). One sees that depending on the accre- 
tion efficiency /, the lifetimes of the PBHs formed during 
the radiation dominated era could exceed the present age 
of the universe ta. 




FIG. 2: The late time evolution of PBH masses (with the same 
initiai mass Mi = 1.5 X lO^^g) for various accretion efficiency 
vafues / = 0, 0.2, 0.4, 0.6. 

On the other hand, integrating the last term in equa- 
tion (1^^ one can obtain the mass of the PBHs which 
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survive beyond the radiation dominated era, given by 



M = 



+ 3o(2„+l)-(^^) ^^^^^ 



(31) 



where Mg = M{te) is obtained by integration of the PBH 
accretion equation ([8|) over the period ti to to be 



1 



Hence, the PBH hfetime is given by 



ievap — i 



1^ 



+ (3a) (2n + l)( 
Further, from equation (|26p . we have 



(32) 



(33) 



= M„ 



M3 



{-©} 

Now using the equations ([55]) and (P5)) . we get 



(34) 



3q! 



2n+ 1 



(^^) M-( 



(35) 



This enables us to invert the PBH hfetime relation in 
order to obtain the formation time for a PBH given its 
time of evaporation t^vapi 



-(a 



(36) 



The above expression is useful for the purpose of evaluat- 
ing the constraints on the initial mass or formation time 
of the PBHs in terms of the observational constraints on 
evaporating black holes at particular eras in this Brans- 
Dicke cosmology. For the present we compute as exam- 
ples the initial masses of the PBHS for two cases: (i) 
PBHs that are evaporating at the present era, and (ii) 
the PBHs that will evaporate when the universe is ten 
times older than its present age. These values are com- 
puted using the analytical result (|36p and displayed in 
the Table I. We have also computed ti and Mi for dif- 
ferent values of the accretion efficiency / from numerical 
integration of equation ([29| . We find that our analyt- 
ical approach gives results that agree up to three dec- 
imal places with the numerical results, thus validating 
the division of the evolution of a PBH into two distinct 
eras dominated by accretion and evaporation dynamics 
respectively, that we have done. 







to 


tevap — 10 X ijQ 


/ 


ti X 10"^^s 


M, X lO^Sg 


ti X 10"^^s 


M, X lO^^g 





2.369 


2.366 


5.105 


5.099 


0.2 


1.658 


1.656 


3.573 


3.569 


0.4 


0.947 


0.946 


2.042 


2.039 


0.6 


0.236 


0.236 


0.510 


0.509 



TABLE I: The formation times and initial masses correspond- 
ing to two specific evaporating eras of the PBHs are displayed 
for several accretion efficiencies. 



IV. CONSTRAINTS ON THE PBH MASS 
FRACTION 

Surviving PBHs at any cosmological era contribute to 
the matter density of the universe at that era. Fur- 
ther, PBHs impact different processes by the end prod- 
ucts of their Hawking evaporation. Various cosmolog- 
ical observations can be used to impose constraints on 
the number density of black holes present during differ- 
ent cosmological eras. These constraints can in turn be 
used for imposing limits of the initial mass spectrum of 
PBHs in various formation mechanisms pertaining to dif- 
ferent cosmological models. In standard cosmology, a va- 
riety of constraints such as coming from considerations 
of overall density, nuleosynthesis, entropy, distortions of 
the CMBR spectrum, and stable relics, have been ob- 
tained [35] . It has been observed that a particular strin- 
gent set of constraints arise from the limits of the 7-ray 
background [111, and also independently from the ob- 
served galactic anti-protons and antideuterons [soj . In 
the present analysis we will just focus on the 7-ray back- 
ground limit in order to obtain bounds on the initial mass 
spectrum of PBHs in BD theory with accretion. 

The fraction of the Universe's mass going into PBHs 
at time t is given by Q 



m = 



PBH 



(37) 



where i^pBH{t) is the present density parameter asso- 
ciated with PBHs forming at time t, z is the redshift 
associated with time t and is the present microwave 
background density having value 10^^. Observations of 
the 7-ray background, as well as those of the antipro- 
tons from galactic sources impose bounds on the present 
cosmological PBH density given by [U, [H, [soj 



PBH 



(t) < 10- 



(38) 



Let us here consider the PBHs that are formed during 
the radiation dominated era, i.e., ti < ti < t^. For them, 
one has 



(1 



(39) 



Using equations ([55)1 and ([55)1 in equation (|57|) . one gets 
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tevap — ^0 


/ 


M^ X lO^^g 


/3(M.) < 





2.366 


5.71 X 10"^® 


1/6 


1.775 


4.95 X 10^2® 


1/3 


1.183 


4.03 X 10-^"^ 


1/2 


0.592 


2.85 X 10"^"^ 


3/5 


0.236 


1.81 X IQ-^"^ 



TABLE II: Upper bounds on the initial mass fraction of PBHs 
that are evaporating today for various accretion efficiencies /. 



a bound on the density fraction given by 



Pit) < 10"^ X - 



(40) 



The above bound pertains to those PBHs which are evap- 
orating at the present era. Again, for ti < t < te, one 

has M = G^^t = Go ^i(l^) • Thus, the fraction of the 

Universe going into PBHs with formation mass Mi is 



/3(M0<10-^x(^ 



(41) 



We can now use the expressions for A/,; (j35p in terms 
of the evaporation time tevap = to obtain the values 
of (3{Mi) corresponding to various values of the accre- 
tion efficiency /.These are displayed in Table IL It was 
shown earlier [38| how the standard constraints on the 
initial mass spectrum are modified in BD theory without 
accretion. Here we observe that increase of accretion effi- 
ciency makes the limits on the initial mass fraction more 
stringent. It may be noted that similar considerations 
would also apply to constraints on the initial mass frac- 
tion [3{Mi) obtained from other physical considerations 
such as those due to entropy or nucleosynthesis bounds. 
The relevant values for the initial PBH masses Mi cor- 
responding to the PBHs evaporating earlier to impact 
entropy production or nucleosynthesis would of course 
be much lower than the values of Mi for which we have 
applied the 7-ray bounds, since these latter PBHs are 
those that are evaporating in the present era. The BD 
dynamics alters the evaporation rate for the PBHs thus 
loosening somewhat the bounds on (3{Mi) as shown by 
Barrow and Carr [IS]. However, inclusion of accretion 
reverses the scenario since accretion is more effective for 
a longer duration for PBHs with smaller which have 
a chance to grow more. As a consequence, ^pBuit) in- 
creases, and the standard constraints on j3{Mi) due to 
nucleosynthesis, entropy, etc., are tightened further in 
the BD scenario with accretion. 



V. CONCLUSIONS 

In this paper we have considered the evolution of 
primordial black holes in Brans-Dicke theory. We use 



the framework of a particular cosmological solution [3j| 
where the gravitational "constant" G can have a much 
larger value in the early universe compared to its present 
strength. We show that accretion of radiation during the 
radiation dominated era is an effective process in this the- 
ory, overriding the mass loss of the PBHs due to Hawk- 
ing evaporation. The cosmological horizon for a PBH 
grows at a faster rate compared to the rate of growth 
of the PBH due to accretion, thus making enough radi- 
ation available for the PBH to accrete during this stage. 
Though the net gain of mass through accretion is not 
significant, it postpones the time for evaporation to take 
over once accretion ceases to be effective at the onset of 
the matter dominated era. The evaporation rate for a 
BD PBH is itself modifed compared to that in standard 
cosmology due the variable G. We show that the life- 
time for PBHs could be enhanced depending upon the 
accretion efficiency /, making the PBHs that are sup- 
posed to be evaporating now (to) live longer by a factor 

of [ 1 ]3/(2n+l) X to. 

The cosmological evolution of PBHs could lead to var- 
ious interesting consequences during different eras. The 
PBHs with smaller masses that are formed during the 
BD field dominated very early era evaporate out com- 
pletely during the radiation dominated era. However, 
those PBHs that are formed later during the radiation 
dominated era survive much longer. These are the ones 
that through their evaporation could impact various cos- 
mological processes such as nucleosynthesis and photon 
decoupling. There exists a variety of observational con- 
straints on the PBHs in standard cosmology |35| . Within 
the context of BD theory Barrow and Carr [24| have eval- 
uated the impact of density bounds on the PBHs evapo- 
rating today on their initial mass spectrum. Here we use 
the observational limits on the 7-ray background [s^ to 
compute the effect of accretion on constraining the ini- 
tial mass fraction of the PBHs. The departure of these 
constraints from those in standard cosmology are quite 
sensitive to the accretion efficiency. 

Finally, it may be noted that there exist other inter- 
esting cosmological solutions for the Brans-Dicke theory 
and its extensions to more general scalar-tensor models 
(40| . In a recent paper one such solution ^ with a time- 
evolving BD coupling parameter w was used to study 
cosmological PBH evolution [25| . It could be worthwhile 
to remap the standard observational constraints on PBHs 
in such models in order to estimate the impact of rapidly 
varying G, as well as possible effective accretion on the 
constraints. Another important issue in PBH evolution 
could be the impact of the back reaction of the PBHs on 
the local background value of the BD field resulting from 
the local change of the density due to the PBHs. Back 
reaction could indeed lead to non-trivial consequences on 
cosmological evolution 4l| , and in the present context it 
might be interesting to see in what way any resulting 
modification could in turn impact the evolution of black 
holes in Brans-Dicke theory. 
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